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Abstract 

The Spinors software is a Mathematica package which implements 2-component spinor calculus as devised by Penrose 
for General Relativity in dimension 3+1. The Spinors software is part of the xAct system, which is a collection of 
Mathematica packages to do tensor analysis by computer. In this paper we give a thorough description of Spinors and 
present practical examples of use. 
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PROGRAM SUMMARY/NEW VERSION PRO- 
GRAM SUMMARY 

Manuscript Title: Spinors: a Mathematica package for doing 

spinor calculus in General Relativity. 

Authors: A. Garcia-Parrado and J.M. Martm-Garcia. 

Program Title: Spinors. 

Journal Reference: 

Catalogue identifier: 

Licensing provisions: 

Programming language: Mathematica. 

Computer: Any computer running Mathematica 7.0 or higher. 

Operating system: Any operating system compatible with 

Mathematica 7.0 or higher. 

RAM: 94Mb in Mathematica 8.0. 

Number of processors used: 1. 

Keywords: General Relativity, Spinor calculus. 

Classification: Relativity and Gravitation. 

External routines/libraries: Mathematica packages xCore, 

xPerm and xTensor which are part of the xAct system. These 

can be obtained at |http://www.xact.es| 

Nature of problem: Manipulation and simplification of spinor 
expressions in General Relativity. 

Solution method: Adaptation of the tensor functionality 
of the xAct system for the specific situation of spinor calculus 
in four dimensional Lorentzian geometry. 



'Corresponding author. 
E-mail address: alfonso@math.uminho.pt 

'Current address: Wolfram Research Inc., 100 Trade Center Dr., 
Champaign IL 61820, USA 



Restrictions: The software only works on 4-dimensional 
Lorentzian space-times with metric of signature 
(1,-1,-1,-1). There is no direct support for Dirac 
spinors. 

Unusual features: Easy rules to transform tensor ex- 
pressions into spinor ones and back. Seamless integration 
of abstract index manipulation of spinor expressions with 
component computations. 

Running time: Under one second to handle and canon- 
icalize standard spinorial expressions with a few dozen 
indices. (These expressions arise naturally in the transforma- 
tion of a spinor expression into a tensor one or vice-versa). 

* Items marked with an asterisk are only required for 
new versions of programs previously published in the CPC 
Program Library. 



1. Introduction 

The concept of spinor plays an important role in 
certain areas of mathematical and theoretical physics. 
Roughly speaking a spinor is a field which transforms 
under a spinor representation of a given symmetry 
group in our system. For example, if we are work- 
ing in a pseudo-Riemannian manifold with a metric of 
signature {p,q) (p represents the number of +1 and q 
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the number of -1 entries in the canonical form of the 
metric), then a natural symmetry group is the group 
which transforms orthonormal frames into orthonormal 
frames. This group is 0(p,q) (SO(p,q) if we restrict 
ourselves to transformations preserving the frame ori- 
entation). The spin group is then the universal cover- 
ing of SO(p,q) which, as is well-known, is Spin{p,q) 
and hence spinors transform under irreducible represen- 
tations of this group. 

The above considerations are completely general and 
they enable us to introduce the notion of spinor field 
in any pseudo-Riemannian manifold admitting a spin 
structure. However, in the case of a 4-dimensional 
Lorentzian manifold (the space-time model in General 
Relativity) a more algebraic approach is desired. This 
approach was pioneered by Penrose lit] where he stu- 
died the main properties of the spinor algebra of those 
spinors arising from the spin group of 5(9(1, 3) and in 
addition he developed a calculus adapted to the parti- 
cular spin vector bundle which one can define in a 4- 
dimensional Lorentzian manifold admitting a spin struc- 
ture. Penrose's spinor calculus revealed very useful in 
certain contexts of General Relativity (GR) where the 
use of tensor methods results in very cumbersome com- 
putations. Perhaps the best known example is the spinor 
formulation of the algebraic classification of the Weyl 
tensor (Petrov classification). The spinor form of the 
Weyl tensor is a totally symmetric 4-rank spinor and it 
is very easy to show that such a spinor can only admit 
four different algebraic types which are in correspon- 
dence with the four distinct Petrov types. 

In this article we describe the Mathematica package 
Spinors which implements the spinor calculus in four 
dimensional Lorentzian geometry as conceived by Pen- 
rose. In this conception, spinors are tensor fields on a 
certain tensor bundle and therefore one can use the gen- 
eral ideas of tensor bundles to work with spinors. In par- 
ticular the notion of spin covariant derivative, the curva- 
ture spinors or the relation between spinors and space- 
time tensors find here a natural formulation. An impor- 
tant part of this formulation is the notion of abstract in- 
dex used to represent tensor fields on any tensor bundle. 
This representation of tensor fields has been adopted in 
the system xAct |01, which Spinors is part of. The sys- 
tem xAct is a system to do tensor analysis by computer 
in Mathematica, both by working with tensors as lin- 
ear combination of basis tensors (component calculus) 
and by working with tensors as symbolic names with 
certain properties like rank or symmetry (abstract cal- 
culus). The system xAct consists of different modules 
tailored for different tasks and Spinors is one of these 
modules. 



Other computer algebra systems support compu- 
tations with spinors. For example in the con- 
text of Particle Physics we may quote the package 
Spinors@Mathematica ^ which can be used in the 
evaluation of scattering amplitudes at tree and loop 
level. The stand-alone package Cadabra handles 
generic abstract spinor quantities in any dimension, with 
emphasis in Field Theory, but no special support for 
General Relativity or component computations. The 
Maple built-in package DifferentialGeometry has exten- 
sive support for component computations of multiple 
types, in particular the NP formalism, but no support 
for abstract tensor computations. Another Maple pack- 
age handling the NP formalism is NPSpinor yl]. 

The paper is organised as follows: in section |2] we 
give a mathematical introduction to spinor calculus. 
The aim of this introduction is to set the notation and 
conventions which are followed by the Spinors imple- 
mentation. Section[3]explains how the Spinors software 
fits into the xAct framework and section |4] presents a 
practical session with Spinors in which the main fea- 
tures of the program are shown by means of practical ex- 
amples. The paper is finished in section |5] where a prac- 
tical computation involving the Nester-Witten spinor 
and the Sparling identity is carried out with Spinors. 

2. Mathematical preliminaries 

In this section we give an overview of the spinor cal- 
culus in General Relativity, following a practical ap- 
proach to introduce the subject and omitting most of the 
proofs (detailed studies can be found in e.g. |0,|3l)- Let 
L be a 4-dimensional real vector space endowed with 
a real scalar product g( , ) of Lorentzian signature and 
let S be a 2-dimensional complex vector space (com- 
plex conjugate of scalars will be denoted by an overbar). 
The vector space S is related to another complex vector 
space S by an anti-linear, involutive transformation. 

The vector space L and its dual L* can be used as the 
starting point to build a tensor algebra in the standard 
fashion. Similarly a tensor algebra is built from S, S 
and their respective duals S*, S . We denote these al- 
gebras by 'X(L), 'X(S) and 'X(S) respectively 0. In this 
work abstract indices will be used throughout to denote 
tensorial quantities: in this way lowercase Latin indices 
a,b, . . . will denote abstract indices on elements of T(L) 



-Strictly speaking only the algebras Xj(L) of tensors r- 
contravariant i-covariaiit can be defined (and the same applies to 
X',(S)). To lessen the notation we will suppress the labels r, s in 
the notation and they will only be made explicit when confusion may 
arise. 
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and capital Latin indices A,B, . . . (resp. primed capital 
Latin indices A',B', . . .) will be used for abstract indices 
of elements in 'X(S) (resp. T(S)). The union of the ten- 
sor algebras 1(S), 1(S) will be referred to as the spin 
algebra and its elements will be called spinors. One can 
also build tensor algebras by taking tensor products of 
elements in T(L), T(S) and T(S). Quantities in these 
tensor algebras will be referred to as mixed quantities 
and they will carry abstract indices of tensor and spinor 
type. All tensor algebras shall be regarded as complex 
vector spaces. 

Since S is 2-dimensional, we deduce that the vector 
space of antisymmetric 2-spinors is 1 -dimensional and 
therefore we can pick up a non-vanishing representative 
eAB which generates such a vector space. We define next 
a spmor e 



^-^B jjjg relation 



^ab£ 



(1) 



where (5^'' is the identity tensor (also known as the Kro- 
necker delta) on the vector space S. Indeed the spinors 
^AB, e^* can be used to relate elements in S and elements 
in S* in the following way 



(2) 



where ^'^ is an arbitrary spinor in S. Hence, the spinors 
and e^^ can be understood as a metric on S (sym- 
plectic metric) and its inverse and the operation shown 
in ^ is the standard "raising and lowering" of indices. 
These operations are extended to the full spinor alge- 
bra without difficulty. In particular we can raise the in- 
dices of £ab getting e^^ = "i^^ and from now on only 
the symbol e will be used for the symplectic metric and 
its inverse. Note also the property 



(3) 



Here the quantity 5^'^ is the Kronecker delta on S and 
is a derived quantity obtained from it by the raising 
and lowering of indices. In particular this implies 6^"^ - 
2. The spinors e^B, e'** and S^^ all have counterparts 
(complex conjugates) defined in the algebra T(S). 

It is possible to relate tensors and spinors by means 
of the soldering form. This is a mixed quantity o-J^' 
fulfilling the algebraic properties 



C AA'^aBB' - ^BA^B'A' 
AA' 



AA'^ b 



cr, 



A A' 



The last of these properties implies that cr^^^, is hermi- 
tian. This is only compatible with the metric signature 
(1,-1,-1,-1). Choosing cr^^, anti-hermitian would 



be only compatible with the signature (-1,1,1,1) |01. 
These properties enable us to relate tensors and spinors 
in the following way 



AiA]...ApA'i, _ a,...fl,, ;,j 
B,B',...B„B„ ~ B,B', 



B„B'^^ai 



AiA' 



(4) 



rj.a,...a„ _ AiA\...ApA' BxB\ ^'iK, a, a,, 

^ L ^ - J o o, o o LT,. ■-■"h " A,A: A„A 



b,...b^ " B,B\...B„B,i 



(5) 



where T ,' " /' is an arbitrary tensor and T „ J,' „ J,' its 

bi...b,i ■' tliB^...B,iB,i 

spinor counterpart 

Another important algebraic property of the soldering 
form is 

O" fl'^A'O" hAC + cr aAC'O'b^A' - gab ^A'C ■ (6) 

This equation is a direct consequence of the irreducible 
decomposition of the product o-g'^A'CrbAC' according to 
theorem[T]below and the algebraic properties of the sol- 
dering form. Starting from (|6]l we can derive formulas 
for the products of soldering forms with all their spinor 
indices contracted (these are useful to translate spinor 
expressions into tensor ones). For example 

A DC' A' 
CraAA'O-b C'Crc CTdD - 

^{iricdab + gabgcd + gcbgda ' gcagdb) , (7) 

where riabcd is the volume form of the metric gab- It 
is possible to generalise this formula for the case of a 
product of more soldering forms. These can be written 
as contracted products of the quantity 

Gabcd = -j(\riahcd + gadgbc - gacgbd + gabgcd)- (8) 

Combining eq. Q and its complex conjugate we obtain 
the spinor counterpart of riahcd, written as follows 

tlabcdCr" AA'Cr^ BB'Cr'^ CC'Cr'' DD' - 
ii^AC^BD^A'D' £B'C' - ^AD^BC^A'C'^B'D')- (9) 

We finish this review about spinor algebra by recalling 
an important result dealing with the decomposition of an 
arbitrary spinor into irreducible parts under the Lorentz 
group M- 

Theorem 1. Any spinor ^Ai -a,,b'^ -b'^, p,q & N can 
be written as the sum of a totally symmetric spinor 
^(A^■■■Ap){B'^■■■B',|) plus terms which are products of the spin 
metric eAB ( or its complex conjugate e^' b' ) times totally 
symmetric spinors of lower rank. 
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2.1. Spinor calculus 

So far all our considerations were algebraic in nature, 
but we can also perform our construction for the case 
of a Lorentzian manifold {M,g) as follows: the con- 
struction performed in previous paragraphs is carried 
out taking as vector space L the tangent space Tp{M) 
of an arbitrary point p e M which is endowed with the 
Lorentzian scalar product g\p. In this way it is possi- 
ble to introduce a complex vector space S,, and a quan- 
tity cr/'^'lp. Now the set S(M) = [jpeAi^p is a vector 
bundle with the manifold M as the base space and the 
group of linear transformations on as the structure 
group. We will call this vector bundle the spin bundle 
and the sections of S{M) are the contravariant rank-1 
spinor fields on M. We can now define the tensor al- 
gebras Tj(7'p(A1)), T^(S/,) and, use them to construct 
vector bundles with Ai as the base manifold. These 
bundles are tensor bundles and we denote each of these 
tensor bundles by 6^'^(A1), where the meaning of the 
labels r, R, s, S is the obvious one. In general we will 
suppress these labels and use just the notation &(M) as 
a generic symbol for these tensor bundles. Sections on 
&{Ai) are written using abstract indices and we follow 
the same conventions as in the case of the vector spaces 
L and S. Sections of any of the bundles ©"f (M) are 
called spinor fields or simply spinors. As usual there is 
a complex conjugate counterpart of this bundle, denoted 
by 6(M). 

Definition 1 (Spin structure). If the quantity crj^ \p 
varies smoothly on the manifold M, then one can de- 
fine a smooth section, denoted by (t^^ . When this is 
the case we call the smooth section crj^ a smooth spin 
structure on the Lorentzian manifold {M, g). 

Clearly a spin structure can be always defined in a 
neighbourhood of any point p € M but further topolog- 
ical restrictions are required if the spin structure is to be 
defined globally (see e.g. Jsl]). 

We turn now to the study of covariant derivatives de- 
fined on the bundles &{M), &{M). Let D,, denote such 
a covariant derivative. Then the operator D„ can act on 
any quantity with tensor indices and/or spinor indices. 
As a result, when Da is restricted to quantities having 
only tensor indices we recover the standard notion of 
covariant derivative acting on tensor fields of M. If Da 
is restricted to quantities having only spinor indices then 
Da is the covariant derivative acting on spinor fields. 
The consequence of this is that the connection coeffi- 
cients and the curvature of D„ will be divided in two 
groups: quantities arising from the tensorial part and 
quantities arising from the spinorial part. The group 



arising from the tensorial part consists of the Christof- 
fel symbols/Ricci rotation coefficients and the Riemann 
tensor of the covariant derivative restricted to the tan- 
gent bundle T{Ai). The group coming from the spino- 
rial part contains the connection components and the 
curvature tensor of the covariant derivative restricted to 
the spin bundle S {M) (or S (At)). We will refer to these 
as the inner connection and the inner curvature respec- 
tively. See Jli 1^ for an in-depth discussion of these 
concepts. 

Definition 2 (Spin covariant derivative). Suppose that 
admits a spin structure (t^^ . We say that a co- 
variant derivative Da defined on &{Ai) is compatible 
with the spin structure cr^ if it fulfils the property 

Dacr,,^^' - 0. (10) 

The covariant derivative Da is then called a spin covari- 
ant derivative with respect to the spin structure cr^^^ . 

Given that any quantity antisymmetric in two spinor 
indices must contain the spin metric as a factor we have 

Da^AB - '^a^AB , '^a = ^^Da^AB , (H) 

where Da is any covariant derivative defined on the bun- 
dle ©(AI). When Da is in addition a spin covariant 
derivative then, combining (fTOl i and ( fTTT i we easily de- 
duce the additional properties 

DaO^AA' = , Dcgab = i'^c + X)gab- (12) 

The last equation implies that any spin convariant 
derivative gives rise to a semi-metric connection ifioll 
when restricted to the space-time tensor bundle. If fur- 
thermore Da has no torsion, then it is known as a Weyl 
connection. The spin covariant derivative gets fixed if 
we demand additional properties on it (see e. g. |[g|) 

Tfieorem 2. There is one and only one torsion-free spin 
covariant derivative Vq with respect to the spin struc- 
ture cr^ which fulfills the property 

^a^AB = 0. (13) 



Acting with such on (|6]l gives 

^agbc = 0, (14) 

which shows that the restriction of V,, to quantities with 
tensorial indices is just the Levi-Civita covariant deriva- 
tive of gab- 
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Consider now any spinor field ^'^ and any spin co- 
variant derivative D,,. Then the commutation of Da, Df, 
acting on is given by t20 

DaDhf - DtDaf - T\,,D,f = F^,Jf (15) 

where T''^^ is the torsion of Da- The mixed quantity 
^oM* inner curvature mentioned above. It is 

antisymmetric in the tensorial indices and it fulfils the 
Bianchi identity 101 

Dv:Fu,\ + VcyA-^- (16) 

The spinor counterpart of the inner curvature is repre- 
sented by Ffjc'DD'A^ ^^'^ decomposed as fol- 
lows 

FcC'DD'AB - XabCD^C'D' + ^ABC'D'^CD- (17) 

The spinors Xabcd and (^abcd' are called curvature 
spinors and they enjoy the symmetries 

XaB(CD) - XabCD , 'i'AB(C'D') - 't>ABC'D'- (18) 

We can also introduce a spinor T^^ g^,^^, representing 
the torsion. Its irreducible decomposition reads 

T^^^ BB'CC = BC^B'C + ^'^^ B'C'^BC, (19) 

where Qaa'bc is the torsion spinor and it fulfills the 
symmetries 

^AA'(BC) — ^AA'BC- (20) 

The inner curvature and the Riemann tensor of D„ are 
indeed related. To find the relation between them one 
computes the Ricci identity for an arbitrary vector V" 
and then particularises it for the special case in which 
y« = cr"^^,^'!^'!'. The result is 

= (K,c%'''' + Fata'''^c'')^f'''c^DD'- (21) 

In the important particular case of a torsion-free con- 
nection which is compatible with the metric (Levi- 
Civita connection) the curvature spinors gain further 
symmetries. These are 

X(AB)CD - XabCD , XabCD - XcDAB , 
'^(AB)C'D' - '^ABC'D'- 

Given these additional symmetries we find that the 
spinor ^abcd' is already in its irreducible form and it 
is called the Ricci spinor. The irreducible decomposi- 
tion of the spinor Xabcd yields 

XaBCD - ^ABCD + M^AD^BC + ^AC^Bd), (22) 



where ^abcd is a totally symmetric spinor called the 
Weyl spinor and A is related to the scalar curvature by 
the formula A = R/24. The curvature spinors and the 
torsion spinor are defined up to a constant scalar factor. 

When working with a spin covariant derivative Da it 
is convenient to introduce the differential operator 

DaA' = (J^AA'Da , (23) 

which enables us to render any expression containing 
spin covariant derivatives as an expression containing 
only spinor indices. Also the commutation Daa'Dbb' - 
Dbb'Daa' can be formally decomposed into irreducible 
parts as follows 

Daa'Dbb' - Dbb'Daa' = ^a'B'Oab + ^abOa'B', (24) 
where 

□as = D(A\A'D^ |g) , Oa'B' = DAiA'D'^B') ' (25) 

are linear differential operators. The action of these op- 
erators on a spinor of any rank is obtained from the 
spinor expression of the Ricci identity of Da and the 
expression of the Riemann tensor in terms of the curva- 
ture spinors. The results for the case of a rank-1 spinor 
are 

□bC^^ = Xd'^BC^" ~ ^'^ '^BC^DA'^^ (26) 

3. The package Spinors and its relation to xAct 

The package Spinors is a Mathematica package 
which implements the spinor calculus as described in 
the previous section. Spinors is part of xAct which 
is a system to do tensor analysis by computer written 
mostly in the Mathematica programming language with 
a smaller part in C. The composition is roughly 16000 
lines of Mathematica code and 2700 lines of C code. 
As of October 2011 the version of xAct is 1.0.3. and 
the complete system is free software available under the 
terms of the GPL license. 

The system xAct is organised as a suite of inter- 
dependent Mathematica packages which can be re- 
garded as different-purpose modules loadable on- 
demand. The packages and the relations among them 
are depicted in figure [1] The arrows indicate which 
packages of the suite a given package relies upon and 
we can see that there are three packages (xCore, xPerm 
and xTensor) which act as kernel for the whole imple- 
mentation. This means that these packages yield the ba- 
sic framework to set up any computation requiring ten- 
sor analysis. In addition the module xperm.c is a piece 
of code in C language devised to speed-up the group 
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The xAct system 



Spinors 



A. Garcia-Parrado and 
J.M. Marlin-Garcia. 



Invar 

'Riemaon tensor Invariants' 

J.M. Martin-Garcia, 
R. Portugal and D. Yllancs. 



SyniManipulator 
"Symmetrized tensor expressions' 
T. BiickdahL 




Harmonics 
"Tensor spherical harmonics" 
D. Brizuela, J.M. Martin-Garci; 
and G. Mena Marugan. 



xPcrt 
"Perturbation theory" 
D. Brizuela, J.M. Martin-Garcia 

and G. Mcna Marugan. 
V- y 



xPrint 
Graphical front-end" 



\Matheniatica Utoisj ^C-languag^: module J 



Finally the formulae relating spinors and tensors had to 
be studied and coded from scratch (this was perhaps one 
of the most time-consuming tasks in the development of 
Spinors). The conclusion of all of this is that spinor the- 
ory is complex enough to develop a new package for the 
xAct suite. 

The package spinors has already been used by a num- 
ber of authors in their research. An example of this is 
the invariant construction of Kerr initial data in |[l5 , 16 ] 
(see also ITtIi for a generalisation of these results). The 
present authors have also used Spinors in the investiga- 
tion of the invariant properties of type D vacuum solu- 
tions of the Einstein's field equations ifisll . 



Figure 1 : Graphical depiction of the xAct packages and their depen- 
dencies. The packages coloured in yellow can be regarded as the sys- 
tem Kernel and have been developed by J.M. Martin-Garcia. 



theoretical computations needed to canonicalize tensor 
expressions lllll . When loading any other package of 
xAct, all the necessary packages are loaded sequentially, 
following the order shown in figure [1] 

The xAct system can be enlarged by adding new pack- 
ages for specific purposes, as long as the dependencies 
just described are kept. Some of the packages shown 
in figure[T]have been already described in the literature. 
This is the case of xPerm 111 ill . Invar lll2lll3ll and xPert 
1 1411 . In addition to these references, each package in- 
cluded in xAct has a documentation file, which explains 
its main features and contains a tutorial. 

The xAct system provides support to define and work 
with general vector bundles on any manifold and dimen- 
sion. This means that the main features of spinor alge- 
bra and spinor calculus are generically built into xAct 
for any dimension. However, in the particular case of 
Lorentzian geometry in dimension four, extra features 
arise as described in section |2] and hence it was nec- 
essary to develop a new package to take care of these 
features. For example the presence of the antisymmet- 
ric metric £ab forces us to be careful with the con- 
ventions to raise and lower indices. Also we need to 
decide which index configuration represents the basic 
Kronecker delta and which index configuration repre- 
sents a derived quantity. In our case all these conven- 
tions were laid in eqs. ©-Ol) and one needs to take 
these particular conventions to the software implemen- 
tation. Another delicate task is the 2-index representa- 
tion of the spin covariant derivatives and their proper- 
ties (eqs. ( |23] )- (l24li). This representation only occurs in 
the spinor calculus and hence we had to code the corre- 
sponding properties explicitly for the Spinors package. 



4. Working with Spinors 

Assuming xAct has been installed one loads Spinors 
in a Mathematica session by typing 

In[l]:= « xAct' Spinors' 

Package xAct ' Spinors ' version 1.0.2, {2011,10,251 
CopyRight (C) 2006-2011, Alfonso Garcia-Parrado Gomez-Lobo 
and Jose M. Martin-Garcia, under the General Public 
License . 



This will load the package Spinors together with the 
other packages of the xAct suite which Spinors relies on. 
These are xCore, xPerm and xTensor (see section [3] and 
figure[TJ. In this work we will only explain the features 
of these packages which are required for our implemen- 
tation and we refer the reader to their documentation for 
further details. 

Next we need to declare a 4-dimensional Lorentzian 
manifold by means of the standard xAct machinery: 

In[2]:= DefManifold[M4, 4, {a, b, c, d, f, h, p)] 
In[3]:= DefMetric[|l,3,0), g[-a,-b], CD] 

The list {a, b, c, d, f, h, pi corresponds to the space- 
time abstract indices which will be used in tensor ex- 
pressions and the list 1 1,3,01 in oef Metric serves to indicate 
the canonical form of the metric tensor g (its canonical 
form contains once +1, three times -1 and zero times 0, 
thus it corresponds to a Lorentzian metric). The sym- 
bol CD represents the Levi-Civita connection compatible 
with the metric g and in addition a number of quantities 
(the Riemann tensor, the Ricci tensor, the Weyl tensor, 
etc) are also created automatically after issuing the com- 
mand Def Metric. 

So far we have used commands belonging to xTensor 
and we have now the set-up necessary to start working 
with Spinors. The first step is the introduction of a spin 
structure. This is achieved as follows 
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In[4]:= DefSpinStructure[g, Spin, (A, B, C, D, F, P|, e, a, 
CDg] 

Several new objects are defined alongside this com- 
mand. These are the spin bundle spin, its abstract in- 
dices |A, B, c, D, F, p|, the spin metric e, the soldering 
form cr and the spin covariant derivative coe compatible 
with both the space-time metric g and the spin metric 
e. The spin bundle spint together with its structures and 
the curvature spinors are automatically defined with this 
command. For example the Weyl and Ricci spinors are 

In[5]:= (PsiCDe[-A, -B, -C, -D] , PhiCDG[-A, -B, -Ct, -Dt]) 
Out [5]= I>P[V]abcd, OmABCtDt-l 

Additional options controlling the displayed form of 
the different quantitites automatically defined can be 
supplied to Def SpinStructure. For example we may add 
the options 

SplnorPref lx-> SP, SpinorMark-> "S". 

The symbol sp will be prepended to the tensor (spinor) 
counterpart of any spinor (tensor) and the string "S" 
will be used in the displayed representation (see be- 
low for explicit examples). From now on it will be un- 
derstood that these options were used in the command 

Def SpinStructure above. 

Primed indices are entered with the "dagger" symbol 
t (entered via the keyboard shortcut ESi3-dg-ESi3). Take 
now the Ricci tensor RiccicD[-a, -u. Its spinor counter- 
part is represented by sp__RicciCD [-a, -At, -B, -Bt] where 

sp was the tag declared through the option spinorPref ix in 
DefSpinstructure. Again, this tag is used to construct the 
symbol defining the spinor counterpart of any tensorial 
quantity, in the way illustrated by this example. In the 
mathematica notebook this is 

In[6]:= SP^RicclCD[-A, -At, -B, -Bt] 
Out [6]= SR[V]AAtBBt 

The linking symbol (entered through the key combi- 
nations EEi3-u[-ESi3) serves to link the symbol chosen to 
represent the spinor prefix with the symbol representing 
the tensor. This linking symbol is stored in the vari- 
able $LinkCharacter which can be freely modified. In ad- 
dition to this, it is possible to modify the default print- 
ing options and obtain outputs similar to the formulae 
described in section |2] This is done as follows for the 
primed indices 

In[7]:= PrintAs[At] "="A"'; PrintAs[Bt] ^ = "B"'; 
PrintAs[Ct] ^ = "C"'; PrintAs[Dt] " = "D"'; PrintAsCFt] ^ = "F"'; 



Also we can modify the printing output of any tensor or 
spinor 

In[8]:= PrintAsEet] "e"; 

The command Decomposition Ca,n be used to find the 
decomposition into irreducible parts of any other curva- 
ture spinor (it is possible to indicate the curvature spinor 
being decomposed as an additional argument). For ex- 
ample 

In[9]:= SP^RieiiiannCD[-A, -At, -B, -Bt,-C, -Ct, -D, -Dt] 
// Decomposition 
□ut [9] = 

^[VltA'B'D'C'eABecD - tlVjABD'C'fDCfA'B' - 2A[ V](eADeBC?A'D' %'C' ^ 
fAC^BD Ca'C'^B'D') ^ 'l*[V]DCA'B'fAB ^D'C + V]ABDCf^A'B' f'C'D' 

The spin covariant derivative cdg can be handled as a 2- 
index covariant derivative ... 

In [10] : = CDe [-F , -Ft] SPsiCDe [-A , -B , -C , -D] 
Out [10]= VFp'i'lVlABCD 

.. or as a covariant derivative in a vector bundle 

In [11] := SeparateSolderingForm [] [%] 

Out [11]= 0-"FF'Va'i'[V]ABCD. 

The command SeparateSolderingForm enables us to trans- 
form spinor indices into tensor ones (cf. 94. II for more 
details about this). 

Spinors are defined by means of Def Spinor, which is 
just a special call to the xAct command Def Tensor (by de- 
fault the option Dagger->compiex is assumcd). 

In[12]:= DefSpinor[K[-A, -At], M4] 

** DefTensor: Defining tensor a-[-A, -At]. 
** DefTensor: Defining tensor A-t[-At, -A]. 

If one wishes to work with Hermitian spinors then this 
is done by using the option Dagger -> Hermitian on the 
command Def Spinor. Under this assumption one has that 
*[-A,-At] is invariant under complex conjugation 

In [13] := K[-k,-kn II Dagger // InputForm 
Out [13]= «[-A,-At] 

The canonicalizer of xAct is Tocanonicai and it can deal 
with the canonicalization of spinor expressions without 
any additional user input. It is beyond the present article 
to explain the workings of the canonicalization proce- 
dure and the reader is referred to |[TT|] and the xTensor 
documentation for additional details. Similarly, the xAct 
command ContractMetric takes care automatically of the 
conventions for raising and lowering indices in spinor 
expressions. We present next some explicit examples 
about these issues 
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In[14]:= Def Spinor [/J [A] , M4] 

** Def Tensor: Defining tensor filkl . 
** DefTensor: Defining tensor /itCAt] . 

In[15] := |/i[-B]e[A,B] [B] e [-B , -A] ,et[Ct,At] 
CDe [-A,-At]a/J[-B]) 

□ut[15]= fBAAl^ e'^'^'^^AA't^B} 

In[16]:= ContractMetricCZ] 
Out [16]= l/jA, /Ja, VA^'/ia) 

In[17]:= I /j[A]/i[-A] , lt[-A, -Bt]//[Bt], CDe[-F, -Ft]aCDe[-A, 
Ft]8PsiCDe[-B, -C, -D, -P]+ CDe[-F, Ft]aCDe[-A, 
-Ft]aPsiCDe[-B, -C, -D,-P] ) 

Dut[17]= l/iAjU*. KAB'/i"'. VfpVa'"'I'[V]bcDP + Vp"" Vaf' VJbCDp) 

In[18]:= ToCanonical [•/.] 
Out [18]= |0, -ka'^'Pb', 0) 

Also the commutation of the spin covariant deriva- 
tives shown in (l24l i is implemented 

In [ 19] ; = CDe [-C , -C t] aCDe [-B , -B t] ®/f [-A , -A t] - 
CDe [-B , -B t] aCDe [-C , -C t] a/f [-A , -A t] 

Out [19] = - Vbb' Vcc'I^AA' + Vcc' Vbb'I^AA' 

In [20] : = SortSpinCovDs ['/. , CDe] 

Out [20] = Ec'B' (□ [ V IcBl^AA' ) + fCB ( □ [ V Jc B-I^AA' ) 

The differential operator \Z\ab can be written in terms 
of the spin 2-index covariant derivatives as illustrated in 
the following example 

In[21]:= BoxCDe[-A,-B]a;f[-C,-Ct] 
Out [21]= □[V]AB/fCC' 

In [22] : = BoxToCovD ['/. , BoxCDe] 
0ut[22]= ^(Vaa'Vb^Vcc + Vba'Va*Vcc') 

In [23] : = BoxToCurvature V/.'/. , BoxCDe] 

Out [23]= - 0|V]ABA'C"fC*' -^[V]dcAB^°C' 

In our previous examples we worked with the spin co- 
variant derivative arising from the Levi-Civita connec- 
tion but we can introduce other arbitrary spin covariant 
derivatives. For example 

In[24]:= Def SpinCovD [nb [-a] , cr, SymbolQf CovD-> |"I","D"), 
Torsion->True] 

As we see in the example the command DefSpinCovD 
shares some similarities with the xTensor command 
DefCovD. In addition to the spin covariant derivative out- 
put symbols, we also need to specify the spin structure 
which the spin covariant derivative is compatible with 
(this is cr in our example). Hence 

In[25]:= |nb [-a] ®o- [b , -A, -At], nb [-a] ecr [-b , A, At] 1 
Out [25]= |(), 0| 



Also, a number of quantities are automatically de- 
fined in addition to the spin covariant derivative nb. 
Since we used the option Torsion- >True the torsion is 
among those and one can work with both its tensor and 
spinor forms. 

In[26] := SeparateSolderingForm [] aTorsionnb [a, -b , -c] 
Out [26] = ST[D]**' BB'CC o^aa' o-b'"'V,<^<^' 

In [27] := PutSolderingFormaDecompositionaZ 

Out [27]= e*Of*'D'(n[D]tDD'B'C'fBC + n[D]D'DBCEB'C') 

In[28]:= ContractMetric9'/. 

Out [28]= n[D]tA*'B'c'eBC +n[D]*'%cfB'c' 

The first step finds the relation between the torsion 
tensor and the torsion spinor and the second step com- 
putes its irreduccible decomposition according to eq. 
(fT9] l. Any spin covariant derivative can be represented 
in single index and two-index notation. 

In[29]:= nb [-A, -A t] 9// [-C] 

Out [29]= ZJaa'Z-'c 

In [30] := SeparateSolderingForm [] [7.,nb] 

Out [30] = o-^AA' i'aA'c 

Finally we remark that it is possible to define a spin 
structure for a metric connection with torsion. In this 
case Def spinstructure defines the torsion spinors automat- 
ically. 

4.1. Relations between tensors and spinors 

One of the strongest points of Spinors is its ability to 
transform tensor expressions into spinor ones and back. 
The transformation rules are ilustrated by ©-(IS]) and to 
work out these expressions in explicit examples we need 
to repeatedly use (O. To illustrate how this works in 
Spinors let us consider the following example: suppose 
that we have the Riemann tensor associated to the Levi- 
Civita connection and we wish to find its spinor form by 
following (|4]i. The procedure is then 

In[31] := PutSolderingForiiiaRiemannCD[-a,-b,-c,-d] 

Out [31]= R[V]„b,j(r'AA'0-''BB"^CC"^'DD' 

In [32] := ContractSolderingForma'/. 
Out [32]= SR[V]aa'bb'cc'dd' 

The Riemann spinor can be transformed back into a 
tensor as follows 

In [33] := SeparateSolderingForm [cr] ['/.] 
Out [33]= R[V];mo-''aa'0-''bb'CT-'cc'CT-''dd' 

In [34] := PutSolderingFormaZ 
Out [34]= R[V]„brf 
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As we see in this simple example a tensor (resp. 
a spinor) is transformed into a spinor (resp. a ten- 
sor) by contracting it with a number of soldering forms 
in the appropriate way. The insertion of soldering 
forms is achieved with the command Put Solder ingForm 
and the elimination of their dummy indices with 

Contract Solder ingForm. The command SeparateSolderingForm 

also inserts a number of soldering forms but unlike 
PutsoideringForm, the cxprcssion on which it acts (a ten- 
sor or a spinor) is automatically replaced by its tensor 
or spinor counterpart. If the tensor or spinor counter- 
part has not been previously defined, then it is created 
automatically. Example: 

In [35]:= Def Tensor [M [-a, -b] , M4] 

** DefTensor; Defining tensor M[-a, -b] . 

In[36] := SeparateSolderingForm [] (SM [-a, -b] 

SpinorDf Tensor :: name : Spinor of M not defined. Prepending 
SP. 

*• DefTensor: Defining tensor SP^M[-q$3,-Qt$3,-Q$5,-Qt$5] . 

** DefTensor: Defining tensor 
SP__Mt[-qt$3,-Q$3,-qt$5,-Cl$5] . 

Out [36]= SMAA'BB'0-.**'o-b'"'' 

All these commands admit a number of options to se- 
lect the tensor (spinor) indices on which one wants to 
act and what tensor (spinor) indices are going to be con- 
tracted. We refer the reader to the on-line documenta- 
tion of each command for the complete list of available 
options. Another related possibility also covered is the 
case in which one has a tensor (resp. spinor) already 
defined in the session and one wishes to introduce its 
spinor (resp. tensor) counterpart. The way in which this 
is achieved is through the command Def SpinorOf Tensor 
(resp. DefTensorOfSpinor). Thcsc Commands allow the 
user to choose the symbols representing the tensor or 
spinor counterparts. For example take the tensor m [-a, -b] 
defined above and suppose that we have not used the 
automatic procedure to define its spinor counterpart de- 
scribed above. Then one can do the following 

In[37]:= Def SpinorDfTensor[SPM[-A,-At,-B,-Bt] , M[-a,-b], 
tr, PrintAs —> "At"] 

** DefTensor: Defining tensor SPM[-A,-At,-B,-Bt] . 
** DefTensor: Defining tensor SPMt[-At,-A,-Bt,-B] . 

The tensor m and the spinor spm are now paired 
to each other For example, if we act on m with 

SeparateSolderingForm [] the SyStCm will USC the SpinOr 

which we have defined above rather than an automatic 
definition. 

In [38] : = SeparateSolderingForm [] (SM [-a, -b] 

Out [38]= AlAA'BB'0-a**'o-b™' 



If the tensor m had had any symmetry, then the symme- 
tries of the spinor spm would have been automatically 
computed. 

When translating spinor into tensor expressions it is 
important to control how products of soldering forms 
transform into tensors. A simple example of such 
a transformation is shown in eq. ^ and products 
of soldering forms with more factors will arise when 
transforming complicated spinor expressions into ten- 
sor ones. The way of computing products of sol- 
dering forms in Spinors is through the command 
ContractSolderingForm. For example, the simplest case is 
the product of two soldering forms 

In[39] := cr[a,-A,-At] o-[b, -B , At] //ContractSolderingForm 
Out [39]= Itr^AB 

The mixed quantity Sct^ab is entered through the key- 
board as sigmao-[a, b, -A, -B] and its square results in the 
tensor Gated introduced in eq. ([8]l. The tensor Gabcd 
will be referred to as the tetra-metric and it is one of 
the quantities automatically defined by Def Metric when 
the manifold dimension is four In this way, if the met- 
ric name symbol is g then Gated is represented by the 
symbol Tetrag and eq. ([8]) by the rule TetraRuie [g] 

In [40]:= Tetrag [-a, -b, -c, -d] 
Out [40]= GgaM 

In [41]:= ■/. /. TetraRuie [g] 

Out [41]= TlEgabrf + ^gailgbc - jgiicgbd + Igabgoi 

The square of Ecr^'^AB is always automatically replaced 
by the tetra-metric 

In[42]:= Sigmao-[-a, -b, -A, -B] Sigmao-[-c, -d, A, B] 
Out [42]= Ggdbac 

The main interest of the tetra-metric is that any con- 
tracted product of soldering forms with no free spinor 
indices can be always expressed as a product of tetra- 
metrics. This is precisely the kind of product which 
arises naturally when translating spinor expressions into 
tensor ones and back. Consider the following example: 
if ^ABCD is the Weyl spinor, we wish to write the scalar 
quantity ^ abcd^'^^^'^ as an expression in terms of the 
Weyl tensor. The Weyl spinor and the Weyl tensor Watcd 
are related through the relation 

"Vabcd = \Watcd a^A^'cr\A'(T'c^' (t'dc , (27) 

and hence the scalar ^ abcd^'^^'^'^ can be computed by 
replacing the Weyl spinor according to dZTl i. Equation 
dZTl i can be written as a xAct rule in the following way 
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In [43]:= WSToWT= IndexRule [PsiCDe [A., B_, C, D J , 
l/4WeylCD[-a, -b, -c, -d]o-[a, A, At](r[b, B, -At] 
(7-[c, C, Ct]o-[d, D, -Ct]] 

Out [43] = HoldPattern [>P[V]^^] : 

Module[|a, At, b, c, Ct, d), \W['7W<lcr'^'^' o-™a'0-''^'^'o-''°c'] 

This construct is called in xAct an index rule. Its diffe- 
rence with a Mathematica (delayed) rule is that dummy 
indices can be included in the right hand side of the 
index rule without caring about the collision of these 
indices with other dummy indices already present in 
the expression in which the replacement is being done. 
Dummy indices will be automatically re-named to avoid 
any index collision. The reader is referred to the xAct 
documentation for further details about this. 

We can use now the rule defined above to find the ten- 
sor expression of any scalar invariant written in terms of 
the Weyl spinor. In our example the actual computation 
runs as follows 



-D] PsiCDe[A, B, C, D] /. 



In [44] := 
PsiCDe [-A, -B, 

Out [44] = 

1 wrv71 wrr7 1 .T-Ji A' bAB' c ^dB C' hCD' ^iD 

In [45]:= ContractSolderingForm [7, , IndicesQf [Spin] ] 
Out [45] = |^g"';?''"'Ci;*,''CgUb» W[ V]crj„,M'[ V]di,i„ 

The option IndicesOf [Spin] USCd in ContractSolderingForm 

indicates that only spinor (dummy) indices in the prod- 
uct of soldering form have to be taken into account in 
the contraction. In this way the final result does not con- 
tain any spinor index and it is thus a tensor expression as 
desired. One can now use the TetraRuie[g] discussed be- 
fore to transform the tetra-metrics into ordinary metrics 
and epsilon symbols (volume elements). 

In [46]:= '/. /. TetraRule [g] ; 

In [47] := ToCaaonicaiaContractMetricSZ 

Dut[47]= iW[V]„b,dW[V]"" + ■,^ifgcdn,W[V]„bn'iy[V]"»='' 



The spinor Haa'bb' is called the Nester-Witten spinor 
and it has the algebraic property "Eaa'bb' + '^bb'aa' - 0. 
Hence, its tensor counterpart, defined by 

- rr rr BB' 
= (^a O-fc "A 

is an antisymmetric tensor and it can be regarded as a 
2-form. 

Theorem 3. The 2-form Tab fulfills the relation { Spar- 
ling identity) 

^V^aTbc] = rjabcf {Z}'^ - \g/^'^ , (30) 

where ^" is the tensor representing the spinor A'^A'^ , Gab 
is the Einstein tensor, rjabcd the volume 4-form (both with 
respect to the space-time metric) and Zabc is a tensor 
fulfilling the "dominant property", namely for any three 
causal future-directed null vectors k", fej, fe^ one has the 
property 

Zabck'{k\k\>Q. (31) 

Proof: We carry out the proof of this result using the 
tools introduced in section |4] (we work in the same 
Spinors session as the one used in that section). First 
of all, we need to define the spinors and tensors inter- 
vening in our problem 

In [48]:= Def Spinor [/I [-A] , M4] 

** Def Tensor: Defining tensor /I [-A] . 
+ * DefTensor: Defining tensor /lt[-At]. 

In[49]:= PrintAsfi^t '^= "I"; 

In[50]:= DefSpinor[H[-A, -At, -B, -Bt] , M4, 
GenSet [-Cycles [|-A, -B), (-At, -Bt)]]] 

** DefTensor: Defining tensor H[-A,-At,-B,-Bt] . 
** DefTensor: Defining tensor Ht[-At,-A,-Bt,-B] . 



In[51] : 



PrintAsiSHt = 



5. Example: The Sparling identity 

As a final exercise with Spinors we show how to use 
the software to derive the Sparling identity. This iden- 
tity has the same information as the Einstein field equa- 
tions and can be formulated either in tensor or spinor 
form. The spinor form of the equation has its origins 
in Witten's proof of the positive mass theorem lfl9ll (see 



also 1 20 1) while the tensor form was found by Sparling 
(see 112 ll l22r for a geometric derivation of this form of 



the equation). 

The set-up is as follows: let A'^ be any rank-1 spinor 
and define the following quantity 



'^AA'BB' = ^('^A'Vbb'/I^ 



'Ib'VaA'^b)- 



(28) 



In [52] : = 

DefTensorOf Spinor [9^[-a, -b] , E[-A, -At, -B, -BT] , cr] 

** DefTensor: Defining tensor ?^[-a,-b] . 
** DefTensor: Defining tensor ytE-aj-b]. 

In[53]:= PrintAsQ^T ^= "T" ; PrintAsOrt ^= "T" ; 

In [54]:= Def Tensor [^ [-a] , M4] 

** DefTensor: Defining tensor ^[-a] . 

We change the default formatting of the volume form 

epsilong 

In[55] := PrintAs [epsilong] ^= "/;"; 

We introduce a short form for the xTensor command 
IndexSolve (see the documentation of xTensor for further 

details about IndexSolve) 
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In [56]:= IndSV [expr_Equal] ; = IndexSolve [expr, FirstOexpr] ; 

We also define the shortcut canonicalization function 
Tc (combination of the xAct commands ContractMetric 

and ToCanonlcal) 

In [57]:= TC[expr_] := ToCanonical [ContractMetric [expr] ] ; 

Also we define a function named EquaiTimes to multi- 
ply by a quantity both sides of an equation and canoni- 
calize the result in just one step. 

In [58]:= EquaiTimes [Equal [lhs_, rhs J , xj : = 
Equal [TC [x Ihs] , TC [x rhs] ] ; 

With all these preparations we introduce the Nester- 
Witten spinor definition, as given by ( l28l l. in our Spinors 
session. 

In [59]:= NesterWittenSpinor = 

E[-A, -At, -B, -Bt] == l/2(-I .it[-Bt] CDe[-A, -At]a/i[-B] + 
I -It [-At] CDe[-B, -Bt]a/i[-A]) 

Out [59]= Eaa'bb' = 5(-Ub' Vaa'-Ib +Ua'Vbb''1a) 

Our aim is to compute the quantity 

In[60]:= AT = CDe [-c] 97^[-a, -b] // Antisymmetrize // TC 
Out [60] = i V.fbc - 5 Vb?\,c + { VcTab 

We transform Tab into the Nester-Witten spinor 
"^AA'BB' and then insert its explicit definition by means 
of the xAct command IndexSolve. 

In[61]:= (.df 1/ SeparateSolderingForm [tr] ) /. 
IndSV [NesterWittenSpinor] ; 

The resulting expression (not shown due to lack of 
space) consists of second covariant spin derivatives of 
Xa and terms formed out of the product aa'Xb, aa~Xb' ■ 
The second covariant derivatives can be eliminated by 
means of the spinor Ricci identity (l24l i and ( l26b . In 
Spinors the procedure for doing this is as follows 

In [62]:= SortSpinCovDs [7. , CDe] ; 

In [63]:= BoxToCurvature ['/. , BoxCDe] ; 

The arising curvature spinors have to be decomposed 
into irreducible parts 

In[64]:= dS = Decomposition[y., Chi, CDe] // TC; 

We split the previous equation into two parts: terms 
containing covariant derivatives of Aa and terms which 
do not contain any covariant derivative. 

In[65]:= (dEl = dE /. CDe [_ ] (8,U_ ] -> 0, dE2 = dE - dEl); 



The aim is now to find the explicit tensor form of each 
part. The first part dBl is an expression which is linear 
in the curvature spinors. We write the curvature spinors 
^aba'B' and ^abcd in terms of the trace-free Ricci ten- 
sor and the Weyl tensor respectively (the rule wstowt was 
defined from eq. dZTl i. see explanations coming after 
that equation.) 

In [66]:= dEl /. WSToWT /. 

IndexRule [PhiCDe [-A., -B., -At., -Bt.], -1/2 TFRicciCD[-a, 
-b] a-[a, -A, -At] o-[b, -B, -Bt]] ; 

Also we need to replace the products A^Aa' by ^a- 

In [67] : = 

•/. /. IndexRule [.1 [B J /lt[Bt.], f[-a](r[a, B, Bt]] // TC 

Finally we eliminate the spinor indices in the previ- 
ous expression. 

In [68] := 

dSl =ContractSolderingForm[7., IndicesOf aSpin] // TC; 

We study next the part containing the covariant 
derivatives of Aa (the expression dS2). This expres- 
sion is a linear combination of spinors of the form 
^bc'^a^cb'Xa' whose tensor counterpart has rank 3. 
This is the tensor we define next. 

In [69]:= Def Tensor [Z [-a, -b, -c] , M4] 

** DefTensor: Defining tensor Z[-a, -b, -c] . 

By definition 

In [70] := PutSolderingForiiiaZ[-a, -b, -c] 

Out [70]= Z;,bi:Cr''AA'0'''BB'0"'cC' 

In [71] := Zrule = 

IndexRule [CDe [-B., -Ct.]ai[-AJ CDe[-C., -BtJ9 /lt[-At.], '/,] 
Out [71]= HoldPattern[Vgf.//l^Vca/^,,/] :^ 
Module[(a, b, c|, Zi,btt'""AA'0"''BB'0"'cc'] 

We use now this rule in the expression for dE2 getting 

In [72] := dE2 = 

ContractSolderingForiii[dE2 /. Zrule, IndicesOf 9Spin] // TC 
Out [72] = 

^iCgbdcfZa'^'^ - ^iC^JbrcdZa**' - ^iG^?adufZb'^^ + ^iCgafL-d Zb"^^ + 
^iGgadbf Zc*^*^ - ^iCgaltdZc'^^ 

We combine now the values just found for dEi and dE2 
and expand the tetra-metrics (see subsection 14. 11 1. The 
final result is 

In [73] := dT == dEl + dE2 /. TetraRuleSg // TC; 

The right hand side of this expression is a compli- 
cated tensor expression of 26 terms. It can be simplified 
though if we compute its double dual (we carry out the 
computation in two steps) 

In[74]:= EquaiTimes [7. , epsilong[-p, a, b, c]] 

Out [74]= /Jp'^'VcT^ib = Z"ap - Z"pa - 35[V]paf" + 3A[V]fp 
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In[75]:= EqualTimes ['/. , epsilongCp, -d, -h, -f]/2] 
□ut [75] = 

Vj?i, - Vr^dh + V|,?df = 

jmavZ"'/ - imthpZ''''„ + 4j;dn,aA[v]f" - \nMpS[v],"e 

In[76];= '/. // TFRicciToRiccl // RicclToEinstein // 
Decomposition // TC 

Out [76] = 

Vd^fh - Vf?-dh + VhTm = j'hiihpZ'/ - jridihpZ"'', - jC[V]aP!;dn,pf" 

This last equation coincides with (l30b and thus we 
conclude its validity. In addition from the spinor ex- 
pression for Zabc we easily deduce the algebraic prop- 
erty ( ISTT i if we express the null vectors ^j, k^^ ^"d k'^ as 
tensor products of spinors of rank-1 (see e.g. theorem 
2.3.6 of Q). □ 

The spinor form of ( l30l l has been used as the start- 
ing point of a proof of the positive mass theorem. The 
rough idea is to prove that the integrals of the Einstein 
and Sparling 3 -form over suitable hypersurfaces extend- 
ing to infinity yield a positive quantity. This is straight- 
forward for the Einstein 3-form if the dominant energy 
condition on the matter is assumed, but it requires more 
efforts for the sparling 3-form. In fact one needs to make 
a special choice of the spinor Aa in order to ensure the 
positivity and there is more than one way of achieving 
this (a good accounnt of the different choices tried can 
be found in Q). 
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